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THE  MICROWAVE  DIAGNOSIS  OF  A  COLUMN  OF  IONIZED  CAS 


J.  R.  Barthel 

ABSTRACT 

* 

The  inverse  scattering  problem  of  deducing  the 
structure  of  a  column  of  ionized  gas  from  its  microwave 
scattering  properties  is  solved  for  a  homogeneous  column 
and  a  column  which  is  sufficiently  dilute  that  the  electrons 
may  be  assumed  to  be  independent  scatters.  Criteria  are 
presented  which  may  be  applied  to  experimental  data  to 
determine  which  of  the  above  models  (if  either)  is 
applicable. 

INTRODUCTION 

A  summary  and  an  extension  of  methods  for  determining  the  signals 
scattered  forward  and  backward  by  a  column  of  ionized  gas,  with  known  properties, 
immersed  in  a  known  incident  plane  microwave  field  is  presented.  A  method  is 
developed  for  solving  the  inverse  problem  of  determining  the  properties  of  the 
gas  column  from  measuring  the  forward  and  backward  scattered  signals  for 
negligible  absorption.  The  direct  and  inverse  problems  are  considered  for 
two  models: 

1,  The  column  of  ionized  gas  is  treated  as  a  homogeneous  dielectric 
constant  is  an  arbitrary  complex  number; 

2.  The  column  is  treated  as  a  large  number  of  individual  electrons 
which  scatter  the  incident  field  independently.  These  electrons 
are  distributed  with  a  number  density  which  may  be  a  strong 
function  of  the  distance  from  the  center  of  the  column  but  a 
weak  function  of  the  distahce  along  the  axis. 


The  theory  Is  applicable  to  the  interpretation  of  measurements  made  of  the 
signal  scattered  forward  and  backward  by  the  ionized  wake  of  a  hypervelocity 
pellet. 

•  * 

SCATTERING  BY  DIELECTRIC  CYLINDERS 

Consider  a  dielectric  cylinder  in  a  vacuum.  The  incident  plane  wave 
propagates  in  a  direction  normal  to  the  axis  of  the  cylinder  as  shown 
below. 


Qforward  RECEIVER 


4 


i 


An  approach  to  this  problem  is  to  expand  the  incident  plane  wave  as  an  infinite 
series  of  cylindrical  waves.1  The  reflected  and  refracted  waves  are  similarly 
represented,  but  the  coefficients  in  these  two  series  remain  undetermined. 

By  applying  the  boundary 'conditions  at  the  surface  of  the  cylinder,  it  is 
possible  to  find  these  coefficients.  However,  for  arbitrary  variations  of 
the  dielectric  constant  K  (|),  an  equation  of  the  following  type2  must  be 
solved  for  each  integer  n  in  the  series  from  -«  to  » 


(|f  T  +  (  k2  K(|)  |  -  n2/;)  f  =0. 
n 


©  Each  choice  of  the  function  K(|)  presents  us  with  a  very  formidable  computing 
task,  even  if  the  ^ries  converges  rapidly.  Thus,  in  practice,  the  study  of 
the  dielectric  cylinder  is  made  only  for  a  very  good  conductor  (refracted 
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O 
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fields  vanish)  and  a  unifora  dielectric  [K(|)  »  constant )j  the  foraer  will 
be  represented  here  as  a  special  case  of  the  latter  by  taking  the  Unit 

We  will  consider  only  the  transverse  magnetic  (TM)  case,  in  which  the  . 
incident  wave  is  polarized  in  the  Zr  direct  ion*  The  transverse  electric  (TE) 
wave,  oolarized  perpendicular  to  the  Z-axls,  has  been  treated  by  Panof sky  • 

and  Phillips1  and  by  Eshlcmn2  as  well  as  by  otters  for  some  of  the  special 
cases  which  we  will  discuss* 

The  general  formula  for  the  electric  field  at  all  points  outside  a  uniform 
dielectric  cylinder  immersed  in  a  uniform,  infinite  plane  C™)  vave  is  the  sum 

*  of  the  incident  and  the  scattered  fields? 

00 

e-lut  |E*l{  [e1*50086  .toSFn(kaiK)]  }  (1) 

* 

where 


J  (Jkka)J*(ka)  -  >/icJ  '  (^Kka)  J  (k*) 
n  n  n  _  n  -  -  f 

aJn  (>/Kka)H^)'  (ka)  -  J*  Jj  (fi.  ka)  (ka) 


J  ,  Hf1)  are  tht  Bessel  and  Hankel  functions  of  the  first  kind,  order  n,  and 
k  =  u/C.  In  all  cases  which  interest  us,  the  point  of  observation  is  far 
enough  from  the  cylinder  that  H^frl)  may  be  represented  by  the  asymptotic 

expression 

•  H^to)  exP  (2) 


which  holds  k£  »  1* 


Here,  Eq.  (1)  becomes 
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It  Is  convenient  to  define  a  backscattering  reflection  coefficient  p^ 
and  a  forward  scattering  coefficient  6uch  that 
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=  jEz|  >f  —  %  * 

0=7 r  '  '  TikH^ 
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E 


0) 


From  Eqs.  (3)  and  (4),  ve  find 

00 

•  t*1'"  [*„(■“»  K>] 


and 


=  T  [Pn^')] 


(5)  # 


(6) 


In  general,  the  evaluation  of  expressions  (5 )  and  (6}  is  very  tedious 
* 

and  must  be  repeated  for  each  given  set  of  values,  of  ka  and  K;  therefore,  it 
will  be  useful  to  examine  certain  limiting  cases* 

Limiting  Values  of  K 

The  dielectric  constant  K  may  take  a  variety  of  forms*  '  For  the  usual 
non-conducting  dielectric  medium,  it  is  a  real  number  >  lj  for  a  conductor , 

.  u,2/^ 

K  =  1  +  —  >  for  an  ionized  gas,  K  =  1  -  .  .  where  w  is  the  plasma 

€u  1  -  iv/,u>  v 

frequency  (=  N.Cele,ctrons/I5  ie._  )  and  v/u>  represents  the  attenuation  by 

m  ec2  • 

e 

collision  and  by  radiation.  t 

In  a  very  good  dielectric  or  a  very  good  conductor,  we  find  |k| 

Then  Eqs.«(5)  and  (6)  approach  the  following  form,  as  |k|  -*•  »  , 


Pb 


-  ■-  i  ^ 


J  (ka) 
n  n 


-00 

00 


Hj^Cka) 


pf  = 


y  X 

L  rrfl 


J  (ka) 


H(x)  (ka) 

_oo  n 


(7) 


(8) 
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(Note:  Eq.  (7)  is  shown  graphically  by  Eshleman,?  pp.  1U  and  112.) 

• 

In  an  •Ionised  gas#near  the  critical  frequency,  |k|«1.  As  !k[  -*  0, 
however,  no  terms  in  Eq.  (1)  may  he  dropped.  It  nay  b<e  seen  by  expanding 
Jn  (s/k  ka)  and  J  '  ('K  ka)  for  small  arguments  that  the  two  teras  in  the 
numerator  of  Eq.  (l)  arid  the  two  terms  in  the  denominator  are  all  of  the 
same  order  in  K. 

Limiting  Values  of  ka  • 

The  number  ka  =  wa/c  is  a  real  number*  We  consider  two  limits: 

1*  fca  «  1#  i.e. ,  the  long* wave  limit* 

2.  Jta  »  If  the  optical  limit. 

For  small  values  of  ka,  we  find  the  following  results! 
large,  so  that  |  «Jk  ka|  <  <  1  as  well  as  ka  «  1,  then 

iir  (K  -  llOca)2 

»b  *  pr  *  r  v  • 

•  » 

3 

If  jx|  »  1  then 


If  |K|  is  not  too 

(9) 

# 


(log  r  r 0.5779;  rz 


fLO) 


Jf  |K|  «  1,  then 

*  -  Pf  -  -  £  <taf  »» 

Hate  that  Eq.  (10)  Is  wild  Only  If  |K|  Is  so  lafge  that  |K|  >af  log 

•  (Ika)  »  1*  © 

• 

For  large  values  of  ka,  it  is  difficult  to  approximate  the  infinite 
series  of  Sq.  (l).  However,  the  laws  of  geometrical  optics  hold  in  this 
limit.  The  reflection  coefficient  from  a  perfectly  reflecting  cylinder 
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*G* 


(for  whirh  jK  |  »  1)  is  readily  found8  to  be 


(12) 


•  • 


ELECTRON  SCATTERING  THEORY 

Consider  on  ionized  column  of  gas  vjth  variable  electron  density* II  (|) 
which  is  immersed  in  a  plane  oierovate  field.  The  incident  vaVe  propagates 
in  a  direction  normal  to  the  axis  of  the  column  os  shown  below* 


W 


* 


The  egression  of  Eshlesaa's*  El*  (64)  for  the  scattered  electric  field  at  * 
the  point  of  observation  for  9  =  0  ttay  be  extended,  for  arbitrary  0  tO  tie 
following  fora:  :• 


a  2f  oo  * 

r  ,r  r*  n(i)  exp  {ik  tfr-scostp-aU-Wsfc^Qt^te  /«% 

K—*r*-  J  J  J  - 9  ;  .  -  -  - — -  j(W) 

me  o  o  »oo  0  •  « 

© 

here  we  have  assumed  that  all  electrons  are  subject  to  an  incident  wave  Of  the 
same  strength.  This  assumption  holds  only  for  sufficiently  dilute  columns 
in  which  v  «  u>.  Note  that  E  and  E  at  the  observer  are  zero  because  of 

*>  x  y  * 

symmetry.  # 
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If  a  back- scattered  signal  is  being  observed  (g  =  0),  and  if  Rq»  a, 
then  the  following  approximate  result  holds 
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"1*0 


2ir 


B  =0 


E*  exp 

2 


, _  a 

5 .  <*>]Nyh  f«u»0v*t)  {ass  a4) 


or 


where 


*  /»UH0( 244) 

*  o 


in 

a 

A  *  -  *-&r*  • 


If  the  fbrvai*d-£tfftttered  signal  is  being  dt  served  J0  =  ir%  and 


Rq  »  a,  then  we  find 


(15) 


=  A 


g=7T 


exp 


_ _  a 

°  o  • 


A_ 

2tt 


exP  [  1  t**o  •?*"*>]  'liiir 


(16J 


where  q  Is  the  number  of  elections  fer  unit  axial  length*  The  total  field 
observed  includes  the  incident  field*  and  my  be  expressed  as 
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exp 


i  (kRQ  -  ut ) 


B=TT 


1+  W~  exp  ('  iT) 

e 


'St  q] 


(17) 


Note  that 


|e  1  < 

E1 

I  z| 

z 

P  =  TT 
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as  expected.  However,  this  effect  Is  not  due  to  absorption,  which  Is  neglected 
but  is  due  to  the  fact  that  the  net  phase  of  the  scattered  signal  opposes  the 
phase  of  the  incident  field  when  observed  in  the  direction  &  » 

Equation  (17)  for  the  field  at  the  forward  observation  point  holds  for 
any  electron  distribution  H(|J»  Equation  (15)  for  the  reflection  coefficient 
*  has  been  evaluated2  for 'several  special  cases:  *  • 


Narrow  Column*  ka«l  ?  *  <i 

k  J,(2ka) 

#  Homogeneous  Cylinder*  N(|)  =  Q, 


o  # 


lea 


Osfbssian  Distribution,  N(|)  =  Nq  exp  (-|2/a2)  : 


(18) 

(19) 


q  exp  t*  (to)2  ]  . 

It  has  also  been  found  that  for  an  exponential  distribution,  N(|)=Ne 


(20) 


-*Va  -s? 


J 


(21) 


1  +  4  k2  a2 


For  arbitrary  functions  N(|),  Eq.  (15)  may  be  evaluated  approximately 
by  expanding  JQ(2k|)  in  series 

a  (2k? )  =  1  -  + 

0  (l!)a  (2  if  (5i)2 


(22) 


a  &  a  a 

\  =  /  »(t)  |d|  -  k* / li|l)l3ds  +£J»(S) sed|  *•••  ■ 


P. 

~K 


This  series  converges  rapidly  if  to  is  not  large. 

If  at  some  part  of  the  column  the  electron  density  is  large  enough  that 
K  <0,  then  the  cylinder  of  radius  |  q  reflects  almost  perfectly  and  may 
^Seated  asVmrtal  cylinder,  if  K  <  0  for  1  <  t0  +  »•  “ere  S  la  the 


8 


skin  depth,  i.e.,  the  depth  Inside  the  surface  Kf  =  0  at  which  the  electric 

field  attenuates  by  a  factor  e_1.  The  region  %  >  |  may  still  be  treated  as 

a  dilute  medium,  but  the  above  analysis  is  not  strictly  applicable  because 

the  field  reflected  free  i  ,  as  well  as  the  incident  field,  acts  on  the 

c 

individual  electron  scatters. 


THE  INVERSE  PROBLEM 

The  problea  of  deteraining  the  structure  of  the  column  from  the  scattered 
signals  Is  more  difficult.  We  wish  to  obtain  the  radius  a  and  the  dielectric 
constant  K,  which  may  be  a  function  of  %;  alternatively,  as  the  problem  is 
posed  in  the  previous  section,  we  wish  to  determine  II (|). 

Only  certain  simple  configurations  may  be  analyzed  successfully.  The 
first  step  is  to  determine  which,  if  any,  of  the  following  simple  models  is 
most  consistent  vith  the  data:  the  dilute  column  of  scattering  electrons,  the 
narrow  uniform  dielectric  cylinder,  or  the  perfectly  reflecting  cylinder. 

Our  choice  of  a  model  must  be  guided  by  a  priori  knowledge  of  the  column,  as 
obtained  from  approximate  calculations  for  the  known  conditions;  this  choice 
must  be  substantiated  by  making  a  comparison  vith  the  experimental  results. 


If  estimates  indicate  that  collisions  are  significant  or  that  q  exceeds 
a  certain  limit  (which  will  be  obtained  later)  then  we  may  postulate  the 
narrow,  uniform,  dielectric  cylinder  model  and  compare  it  with  experiment. 
From  Eq.  (9)  ve  obtain 


lira2  Wp2 

^  =  Pf  =  c2(l  -  iv/u) 


(257 


We  may  determine  (a  up)  and  v  by  applying  Eq.  (25)  at  two .frequencies.* 

The  model  becomes  convincing  only  if  there  are  a  number  of  measuring 
frequencies,  so  that  the  result  may  be  checked  several  different  ways. 

Note,  however,  that  we  may  apply  Eq.  (25)  only  at  frequencies  for  which  ka 
is  somewhat  less  than  unity.  The  indeterminacy  may  be  resolved  if  a  frequency 
can  be  found  which  is  low  enough  that  the  skin  depth  &  is  much  less  than 


*Note  that  (a  u)p)2  =  q  x  known  constant,  i.e.,  we  then  know  the 
number  of  electrons  per  unit  axial  length. 


the  column  radius}  we  find  5^^  —  c/Wp,(as  u-*0).  If  c/wp  «  a,  we  may 
use  Eq.  (7)  to  find  a,  since  the  column  is  then  effectively  a  perfectly 
reflecting  cylinder}  if  not,  the  indeterminacy  will  remain  unless  some 
independent  way  of  finding  a  or  up  is  found. 

If  the  column  is  not  sufficiently  uniform  for  the  above  method  to 
succeed,  then  it  appears  that  a  simple  microwave  analysis  of  the  column 
structure  Js  not  feasible  unless  the  dilute  column  model  (independent  * 
scattering  electrons)  may  be  applied.  However,  this  model  applies  only  when 
absorption  (by  collisions  with  neutral  particles  or  by  electron  resonance 
(u  <  u  ))  is  small.  For  this  model  Eq..  (17)  leads  to  the  following  relation 

for  the  forward- scattered  power  ratio* 


i.i.i  fef 

P._  2  Vm. 


where  P  refers  to  the  absence  of  the  ionized  column.  The  solutions  of 
Eq.  (2lf)°for  q,  from  the  measured  values  of  Pf/P^*  should  give  nearly  the 
same  answer  for  each  k  >  up(max)/c,  if  we  are  t*be  Justified  in  using 
the  dilute  column  model. 

•  If  application  of  Eq.  (24)  for  various  k  Justifies  the  dilute  column 

model,  we  may  then  proceed  to  determine  N(|)  by  solving  Equation  (15) 

00 

•  Pb(k)=A  / N(|)  J0(2k|)  |d|.  (2 

*  o 

The  solution  is  found  by  applying  the  Fourier-Bessel  integral3  which  states 


=  fjm  (kz)  kdk  J  f(|)  Jm(kS)Sd|  (m  >  -  JO 


Note  that  if  q  exceeds  a  certain  value,  which  need  not  correspond 
to  w  >*<*),  there  is  a  focusing  effect,  i.e., 

p 
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After  applying  Eq.  (26),  we  find 


J  (2k» )  (2k)  2dk 
o 


N(x)  J  (2kx)  xdx 
o 


« 

J0(2kt)  Oh  (It)  kdk. 


(2T) 


♦ 


There  is  a  question  of  the  validity  of  the  application  of  Bq.  (27). 

There  is  a  frequency  u',  or  a  wave  number  k ",  below  which  the  Integrand 

does  not  have  the  form  which  corresponds  to  a  dilute  column  (u' Is  the  plasma 

frequency  corresponding  to  N^).  If  the  region  k  <  k'  does  not  make  Ihe 

dominant  contribution  to  the  integral,  we  may  resolve  this  question  by 

employing  the  following  scheme.  First  we  determine  a  trial  function  N(|) 

by  Eq.  (27)  by  using  the  experimental  values  of  f^(k).  ve  consider 

the  vfclue  N  >  which  corresponds  to 

max  • 


b) 


*2  _ 


u^(max)  = 

P 


e2li. 


max 


m  e 


(k'c)2 


(28) 


We  assume  N  =  K(0),  for  simplicity  in  developing  the  scheme,  although 
max 

N  4  N(0)  is  treated  by  .a  straightforward  extension  of  the  following, 
max 

Now  notice  that 

00 

N(0)=XJ  f^(k)kdk  =  Nmax 

o  • 


(29) 


We  wish  to  adjust  p^k),  for  k  <k',  to  the  form  it  would  have  if  the  column 
continued  to  act  as  a  number  of  independent  scattering  electrons,  i.e.,  we 
wish  to  remove. the  effect  of  absorption  and  reflection  from  regions  in  which 
to  >  u.  We  know  N  ,  we  find  k ',  and  then  we  analytically  extend  the 

■n  JU&X  *  ,  . 

function  ^  (k)  (for  k  >  k')  to  the  interval  0  <  k  <  k  .  In  so  doing,  we 
may  utilize  the  fact  that  for  the  dilute  column  model  with  ka«  1,  pb  =  qA/27T, 
where  q  is  estimated  by  using  the  trial  function  N(|)  in  the  expression 

00 

q  =  27 T  [  N(|)|d|  . 
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By  using  the  adjusted  function  ^(k),  «e  then  repeat  the  above  procedure, 
and  again  find  a  trial  function  H(|),  find  k',  and  adjust  ^(k)  for 
k  <  k'.  We  Bay  continue  to  repeat  the  procedure  until  convergence  is 
evident.  However,  it  is  doubtful  that  the  scheme  will  lead  to  accurate 
values  of  H(fr)  unless  the  function  to  which  H(|)  converges  is  not  greatly 
different  from  the  result  of  the  first  correction}  i.e.,  we  oust  restrict 
k'  to  be  small  enough  so  that  the  major  contribution  to  the  Integral  for 
H(|)  at  each  |  is  In  the  region  k  >k'. 

We  may  estimate  the  severity  of  the  above  restriction  by  considering 
the  case  *'  a  «  X;  than  ^  <  *')  =  iV&r  t*.  P»H-  *b0Ve 
restriction  requires  that 


N(0)  » 


u 


k' 


Pb 


k  d  k 


_  Sa  fk')2  qe2N(0) 


me  c 
e 


00) 


or  q  <  <  1014  elect  rons/me  te^  where  ve  hate  employed  Eqs.  (28)  and  (29). 
Therefore,  ve  may  expect  successful  appUeatlon  of  Eq.  (27)  to  the  inverse 
problem  only  If  approximate  calculations  for  the  prevailing  condition  in 
the  column  indicate  that  q  is  considerably  leas  than  10»  electrons/meter 
of  column  length.  This  restriction  may  he  unnecessarily  severe,  because 
the  value  of  a  assumed  for  k  <  k'  is  the  maximum  that  any  of  the  models 
in  Eqs.  (18)  to  (21)  may  attain,  nevertheless,  if  ve  assume,  for  example, 
a  Gaussian  distribution,  the  above  restriction  leads  to 

e'(to)2  k  d  k  =N(0)[l-e"(k  a)2] 


N(0) 


k  d  k  = 


=  i+a 


2ir 


or  (kra)2  «  1 


which  vill  lead  to  the  same  restriction  as  Eq.  00). 

There  is  less  restriction  on  the  validity  of  the  scattering  electron 
theory  for  the  direct  problem.  If  the  plasma  frequency  wp  in  I*88  than  “ 
everywhere  in  the  column,  and  if  »/«  « 1,  i8  valid  for  the 

direct  problem  no  matter  how  large  q  may  be.  Eshleman  has  shown  that  for 
a  column  with  mp  >  u  and  kro  <  1,  where  rQ  is  defined  by  q  -  T 
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then  it  is  required  that  q  «  XO14  electrons/raeter  to  validate  the  theory; 
if  (krQ)  >  1  and  w  >  u,  It  la  required  that  q  «  (krQ)£  1014  electrons/meter. 
The  core  stringent  requirements  for  validity  of  the  theory  in  the  inverse 
problea  are  caused  by  the  necessary  consideration  of  very  small  as  veil  as 
very  large  values  of  «. 


SUMMARY  AND  CONCUUSIOKS 

We  have  first  summarised  results  derived  elsewhere  for  microwave 
scattering  free  collects  whicli  are  cither 

1.  uniform  dielectrics  (K  cay  be  complex),  or 

2.  sufficiently  dilute  that  the  electrons  scatter  independently. 

These  results  were  titan  extended,  for  the  second  case,  to  forward-scattering. 
It  was  found  that  the  ratio  of  the  power  of  the  forward- scattered  signal 
to  its  value  in  the  absence  of  a  scattering  column  may  exceed  unity,  i.e., 
there  is  a  focusing  effect;  it  wus  also  found  that,  neglecting  collisions, 
this  power  ratio  depends  only  on  the  number  of  electrons  per  unit  length 
of  the  column. 

The  inverse  problem  of  determining  the  column  structure  from  its 
scattering  properties  was  also  studied,  for  both  cases.  In  Case  (l.)  we 
may  determine  at  the  least,  the  collision  frequency  v  of  the  electrons 
and  the  number  of  electrons  q  per  unit  axial  length;  if  the  column  reflects 
perfectly  at  sufficiently  low  frequencies,  we  may  measure  the  column 
radius  a,  and  hence  deduce  the  electron  density  N  from  the  value  of  q. 

In  Case  (2.)  it  was  show  that  the  electron  density  N(|),  which  may 
vary  radially,  may  be  determined  from  the  reflection  coefficient  by 
applying  the  Fourier-Bessel  integral,  provided  q  is  considerably  less 
than  1014  electrons/meter. 
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